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We investigate theoretically the linear and nonlinear conductance through a nanostructure with two-fold 
degenerate single levels, corresponding to the transport through nanostructures such as a carbon nanotube, 
or double dot systems with capacitive interaction. It is shown that the presence of the interaction asymmetry 
between orbits/dots affects significantly the profile of the linear conductance at finite temperature, and, of 
the nonlinear conductance, particularly around half-filling, where the two-particle Kondo effect occurs. 
Within the range of experimentally feasible parameters, the S U(4) universal behavior is suggested, and 
comparison with relevant experiments is made. 
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1. Introduction 

Electrical transport through nanostructure devices or a 
quantum dot is fundamentally influenced by the presence of 
the Coulomb interaction on the dot, where along with the 
Coulomb blockade phenomena, the Kondo effect, a paradigm 
of the strongly correlated phenomena, has been observed. '' 
The great advantage of these nanostructure systems is the 
ability to control various relevant parameters to regulate 
many-body effect. Toward future applications of quantum 
electron devices, interest in a role of many-body effect of 
the transport phenomena has been rekindled and studies have 
been made actively. 

In the absence of magnetic field, the levels of a quantum 
dot are spin degenerate; the Kondo effect gives rise to the con- 
ductance enhancement at low-temperatures in an odd electron 
number on the dot, as a result of the interaction between a lo- 
cal magnetic moment of the dot and conduction electrons. The 
spin degree can be replaced by another degenerate degrees of 
freedom. Indeed, two quantum dots with capacitive interac- 
tion can be labeled by an index i = 1,2, which can be taken 
as another realization of the Anderson impurity model; 2 ' the 
orbital Kondo effect appears. 3 ~ 10 ' Other notable examples are 
vertical quantum dots, 11 ' or carbon nanotube dots. 12-15 ' All 
of these systems accommodate one, two, three, four electrons 
in the topmost shell, and the Coulomb blockade peaks show 
a clear sign of the fourfold pattern. When the spin and or- 
bital degeneracies are simultaneously present, the further en- 
tanglement between spin and orbital degrees of freedom oc- 
curs, giving rise to the S U (4) Kondo effect with an Kondo 
temperature one order of magnitude higher than the standard 



case 



3,4,9,10,14-18) 



Though most of studies of the systems with doubly de- 
generate orbitals have been focusing on odd-number valleys 
(Nd ~ 1 or 3), it has been recognized recently that a simi- 
lar Kondo-like enhancement (two-particle S U (4) Kondo ef- 
fect) can occur on even-number valley (Nd ~ 2) as well 19 ' 20 ' 
(See also ref. 21 for the singlet-triplet Kondo effect.) How- 
ever, such enhancement at even-number valleys does not seem 
ubiquitous, observed in some experiments but not in oth- 
ers. The conductance of carbon nanotube dots fits well with 
G oc sin 2 (7rA^/4) at very low temperature, 15 ' which shows the 
Kondo enhancement both at odd and even valleys; but in two 



coupled dots, the honeycomb-like structure as a function of 
the two gate voltages shows the enhancement at Nd = 1 but 
not at N d = 2. 

In this paper, we will examine the quantum transport 
through a dot with doubly degenerate orbitals with different 
interaction strengths within the orbits and between orbits. We 
will show that such interaction asymmetry, if small, is im- 
portant to understand the profile of the linear conductance at 
finite temperature, and, of the nonlinear conductance. Resort- 
ing to the universality of the topmost shell, the S U (4) Ander- 
son model, we analyze how the conductance gets enhanced or 
suppressed at finite temperature or by finite source-drain volt- 
age. Though the conductance profiles observed in several ex- 
periments may seem nonuniversal at first glance, we will find 
that we can explain them systematically by the SU(4) An- 
derson model universality, once taking proper account of the 
presence of the interaction asymmetry between orbits/dots. 

We are concerned particularly with the universal depen- 
dence of the temperature and/or the bias voltage of the 
conductance. Naively, one may expect the singlet-triplet 
Kondo effect at half-filling when the interaction asymmetry is 
present. On the other hand, even with the interaction asymme- 
try, the renormalization group (RG) flow is known to drive the 
system toward the S f/(4)-symmetric strong coupling point, 
either at quarter-filling or at half-filling. 19 ' 20 ' Because of it, 
it is unclear what is the relevant symmetry of the system 
once interaction asymmetry is present; the unitarity limit of 
the conductance enhancement at quarter-filling cannot distin- 
guish between them. In the standard spin-degenerate S U (2) 
Kondo systems, experiments have confirmed universal scal- 
ing regarding temperature 22 ' as well as the bias voltage at odd 
valleys. 23 ' 24 ' Since the absence of the exact SU(4) symmetry 
exhibits the behavior at finite temperature and/or at finite bias 
voltage, we will try to identify the relevant universality, S U(4) 
or S U(2) x S U(2), by examining the universal scaling of the 
linear/nonlinear conductance. We will show numerically no 
crossover observed in the conductance for a experimentally 
relevant parameters U'/U = 0.4 ~ 1.0 and U = 10 ~ 20T; the 
Kondo-like enhancement is controlled by the universal behav- 
iors of the S U (4) Anderson model only with the characteristic 
temperature renormalized. 
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Outline of the paper The paper is organized as follows. In 
Sec. 2, we introduce the model of quantum transport through 
a nanostructure with two-fold degenerate orbits, and main ap- 
paratus of our approach, finite interaction slave-boson mean 
field theory; a brief theory of the linear/nonlinear conductance 
is summarized. In Sec. 3, we present our theoretical results 
both for the linear and nonlinear conductance. We see the 
importance of the interaction asymmetry on the conductance 
profile. A direct comparison with numerical RG data is also 
made. After introducing the characteristic temperature T* in 
Sec. 4, we discuss our theoretical results in light of T*, and 
universal scalings of linear and nonlinear conductance is ex- 
amined in Sec. 5. After discussing an issue of T* renormaliza- 
tion versus a crossover in Sec. 6, we compare and discuss our 
results with typical experiments where the interaction asym- 
metry is expected to play an important role in Sec. 7. Finally 
we conclude in Sec. 8. Part of the results on the linear con- 
ductance has been reported earlier. 25 ' 

2. Model and Formulation 

2.1 Quantum dots with doubly degenerate orbitals 

Relatively large energy scale of the Coulomb interaction 
allows us to focus on the topmost electron shell. Our ba- 
sic assumption is that within the universality of the topmost 
shell, the Anderson impurity model can describe appropri- 
ately quantum transport phenomena for any value of N d . We 
model the topmost shell of the dot as degenerate in orbits 
i = 1,2. Such degeneracy may come from the coupled-dot 
or a unique electronic structure of carbon nanotubes. 

On the dot, an electron interacts with an electron in the 
same orbit i — 1, 2 by U or in the different orbit by U'. The 
total Hamiltonian is given by 

H = H D + H L + H T , (1) 

where the dot H D , the noninteracting leads H L , and the cou- 
pling between leads (a = L, R) and the dot H T are defined 
by 

H D = ^ Ed h kT + ^ t/%% + U'hih 2 , (2) 

icr i 

H i = ZZ fe " ^a) C^Cctor, (3) 
a kicr 

*r = 2Z( ftc LA +h - c -)- ( 4 ) 

a kicr 

Here the number operator of the dot is defined by n, = 
Ho- "<<x = Za- dl^div and the gate voltage Sd control the av- 
erage electron number on the dot Nd = «;) continuously 
from to 4. In the calculation below, we assume constant den- 
sity of states p a of the lead a within the energy band (-D, D), 
and use r = 2Z a ^a = 2Za n Pa\tk\ 2 as a coupling parameter 
between the leads and the dot. To investigate the nonlinear 
conductance in the presence of finite source-drain voltage V, 
we choose the chemical potentials [i a as fi L = -fi R = eV/2. 

In the case of U' = U, the total Hamiltonian H exhibits the 
full S £7(4) symmetry. The states (n u n 2 ) = (1,0) and (0, 1) 
are four-fold degenerate at quarter-filling (Nd = 1); n 2 ) = 
(2, 0), (1, 1) and (0, 2), six-fold degenerate at half-filling (Nd = 
2). 

When one breaks the S U(4) symmetry by decreasing U', 
the four-fold degeneracy at quarter-filling is unbroken, but 



the six-fold degeneracy at half-filling is broken. The effect of 
the interaction asymmetry appears as shifts of the Coulomb 
blockade peaks, those ats d = 0, -U', -U-U', and -U -2U'. 
However, the effect of the interaction asymmetry does not 
stop here; it gives a substantial effect in the Kondo effect par- 
ticularly at half-filling. This simple argument also indicates 
that the Si/ (4) symmetry at half-filling is more vulnerable 
than that at quarter-filling. We will show below that this is 
indeed the case. 

2.2 Finite -interaction slave-boson mean-field theory 

We analyze the model by an extension of the Kotliar- 
Ruckenstein formulation of slave-boson mean field theory 
(KR-SBMT), 26 ' where a bosonic field is attached to each type 
of local excitations. 27 ' The approach has several advantages 
that other slave-boson cousins miss. By retaining finite in- 
teraction, it can treat the mixed-valence regime; it applies to 
nonlinear transport through a quantum dot; 28 ' it reproduces 
Fermi liquid behavior at T = with satisfying the Friedel 
sum rule; it can evaluate the linear/nonlinear conductance for 
the entire range of the gate voltage systematically. The KR 
formulation of SBMT is believed to be an analytical yet reli- 
able non-perturbative approximation up to the Kondo temper- 
ature, agreeing successfully with numerical RG methods and 
experiments 28 ' 29 ' (see also Sec. 3). 

When the KR-SBMT is extended to the dot with doubly 
degenerate orbitals, 16 bose fields are needed associated to 
each state of the dot: e for the empty, p ia for one electron with 
orbit i and spin cr, x, for two electrons on the same orbit i, y sm 
for two electrons at different orbits with total spin (s, m), h i(T 
for three electrons with a hole on icr, and b for fully occupied 
state. Below we briefly summarize our approximation scheme 
(see also ref. 30). 

As a standard treatment of the Kotliar-Ruckenstein formu- 
lation, we can rewrite the Hamiltonian eq. (1) in terms of these 
newly defined fields by adding the Lagrange multiplier terms 
H A = A m (I - 1) + AfJ [d] a d icT - Q io ) to respect two con- 
straints: the completeness condition, 

I = e 1 e + p'Lpitr + J] x]xi 

icr i 

+ 2 yLy»* + J] h ^ hi& + bfb = 1 > ( 5 ) 

sm icr 

and the number correspondence between fermions and newly 
defined bosons, 

QiT = p'lpia + x]xi + y\ M yi,2<r + ^(yl Q yoo + y\ y lo) 
+ h]Ji & + JXAr< +tfb = dld Ur . (6) 

cr' 

Here and hereafter, we adopt the notation 1 = 2, 2 = 1, f =1, 
I =f. As a result, we can rewrite eq. (1) exactly as 

H = H L + J^s d d'ld io - + Uj^ xjxi + U'Yy\ m y sm 

icr i sm 

+ (U + 2U') 2 h\ a h i(T + (2U + AU')tfb 

icr 

kicr 
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where z io - = (1 - Qio-) * [e l Pia- + p^Xi + P^yiM + P^Jym + age V by 

y w )/2 + xihw + CyJo + y\ )kfrll + y[ 2a .h 1(r + h\b\ QTJ denotes 
a renormalization of one-particle annihilation. 

By applying the mean field approximation, we replace all 
the boson fields by their expectation values. We determine 
those expectation values self-consistently by solving equa- 
tions of motion and constraints at each gate voltage e d , the 
bias condition V and the temperature T. By adopting a sym- 
metric solution regarding the orbit ;' and the spin cr, we 
should solve the following self-consistent equations (omitting 
spin/orbit indices), 



4 d ^l M + A^e = 0, 
de 

d\nz 
Sinz 

dx 

dlnz 



M + (A m - A (2) )p = 0, 
M + (U + A m - 2A (2) )x = 0, 



8y 
dlnz 



(8) 
(9) 
(10) 
(11) 



dh 
,dlnz 



M + (U' + A (1) - 2A (2) )y = 0, 
M + (U + 2U' + A m - 3A (2) )h = 0, (12) 



M + (2U + 4U' + A m - AA {2) )b = 0, (13) 



as well as two constraints 

e 2 + 4p 2 + 2x 2 + 4y 2 + 4h 2 + b 2 = 1, (14) 

N d /4 = p 2 + x 2 + 2y 2 + 3h 2 + b 2 . (15) 

Here we have introduced (non-equilibrium) Green functions 
Gj(t,f) = i{d] a .{f)d i(T {t)), its Fourier transformation G^(to), 
a ndM = J^(co-s d )G<(co). 

2.3 Current formula and nonlinear conductance 

On determining these auxiliary parameters self- 
consistently at each temperature and each gate voltage, 
the system reduces to the renormalized resonant level model 

#eff = H L + y| Sd hi + {t k c\ kkT d icT + h.c. ), (16) 

i akio- 

with the effective dot level Sd — + /l (2) and the effective 
hopping tk = ztk, which corresponds to the effective peak 
width f = Fl + Tr = zZa n Pa\tk\ 2 - The effective Hamilto- 
nian conforms to Fermi liquid description at low temperature, 
the strong-coupling fixed point. Note that interaction effect 
is considered only through Sd, h and we ignore a remaining 
(renormalized) interaction term between quasiparticles. 

The form of H e <g enables us to find the linear/nonlinear con- 
ductance by the Meir-Wingreen formula: 31 ' 



/ = 



16ef L f fi 



[/ t M-/«H 

(co - s d ) 2 + f 2 



(17) 



4G r 



Im^(i+^ s )-^(i+^ L )], (18) 



where fdai) = l/[e (t0 fla ^ T + 1] is the Fermi distribution func- 
tion of the lead, G = (4Y L Y R /Y 2 )(e 2 /h), £, = (e d - p a - 
iY)l{2nT), and if/ is the digamma function. By this, we can 
readily evaluate the nonlinear conductance at finite bias volt- 



0(T,V)=§. 



(19) 



We should bear in mind, however, that the renormalized pa- 
rameters s d and f still depend on V as well as T and interac- 
tion strengths. 

At T = 0, the electron number on the dot is equal to Nd = 
2 - £ arctan(ej/f). This immediately shows that the linear 
conductance at zero temperature G(0, 0) satisfies 

G(0, 0) = 4G sin 2 (7r^/4). (20) 

This is the Friedel sum rule of the S U(4) Anderson model. 
The formula immediately indicates that the zero-bias con- 
ductance in the unitarity limit approaches G(0, 0) — > 2Go at 
N d = 1,3, but G(0, 0) -> 4G at N d = 2. The observation of 
the above dependence of G(0, 0) on continuous Nd is a hall- 
mark of the S U(4) Anderson model; a recent experiment 15 ' 
and calculations by numerical RG 32) confirmed it beautifully. 

3. Linear and nonlinear conductance 

In this section, we present our main numerical results of 
linear/nonlinear conductance as a function of the temperature 
and/or bias voltage. The Friedel sum rule indicates the zero- 
bias conductance at T — does not depend on interaction 
asymmetry, taking a universal form, eq. (20). Accordingly, 
the absence of the exact S U(4) symmetry appears at finite 
temperature and/or bias-voltage; we will find that the effect is 
large enough to modify the conductance profile substantially, 
providing a characteristic 'dip structure' around half-filling. 
For all the calculations, we fix a parameter U/F — 20 as a 
typical value of strong correlation in experiments, and we take 
U as a unit energy scale, if needed. 

3.1 G(T, 0) of symmetric interaction U' = U; comparison 
with NRG results 

As a first step, we compare the zero-bias conductance 
G(T, 0) obtained by our KR-SBMT approach with the NRG 
calculations available for the SU(4) Anderson model with 
symmetric interaction U' = U. 3T> We find that our results re- 
produce NRG data quite well. 

Figure 1 (a) demonstrates a distinctive feature of the con- 
ductance profile in large U /F region, the Kondo enhancement 
both at Nd ~ 2 and at N d ~ 1 ; with increasing temperature, 
the four-fold Coulomb blockade peaks appear. Regarding the 
mixed-valence regime with a smaller value of U/F, four-fold 
peaks merge to form one big peak showing the Kondo en- 
hancement (see Fig. 1 (b)). As was claimed already, 32 ' these 
temperature evolutions by S U(4) symmetric Anderson model 
agree very well with what is observed either at V g ~ 3.9V or 
at V g ~ 5.3V in an experiment. 15 ' 

Our KR-SBMT results agree even quantitatively well with 
NRG results, but with one catch: In comparing with NRG 
data, we need to resort to a heuristic prescription choosing 
the twice as large value of U/F for the SBMT results. 33 ' In- 
deed, the same problem has been prevailing in the single im- 
purity Anderson model; we can attain quantitative agreement 
between effective theories (the KR-SBMT 29 ' or the functional 
RG 34 ') and the NRG method only if we choose the twice as 
large value of U/F for the former. Figure 1 shows compar- 
ison by using this heuristic prescription, and they agree very 
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Fig. 1. Comparison between the KR-SBMT and numerical RG results. 32 ' 
Temperature evolution of linear conductance G(T, 0) is plotted as a func- 
tion of the gate voltage for the fully symmetric SU(4) Anderson model, 
(a) SBMT results (lines) of U/T = 20 are compared with NRG data (sym- 
bols) of U/T = 10; (b) SBMT results (lines) of U/T = 10, with NRG data 
(symbols) of U/T = 5 (see the text for details). 



well in low temperature, roughly up to the Kondo temperature 
scale. 

3.2 Effect of interaction asymmetry on G(T, 0) 

Figure 2 demonstrates typical temperature evolutions of the 
conductance profile as a function of the gate voltage with 
varying asymmetry parameters (a) U' /U - 1, (b) 0.8, and 
(c) 0.6. Overall widths of conductance profiles are deter- 
mined by U + 2U'. We notice, even with a small asymmetry 
U'/U = 0.8, a characteristic dip structure of G(T,0) devel- 
oping clearly at half-filling (Nd = 2) at finite temperature, in 
contrast with a small deformation around quarter and three- 
quarter-fillings (Nd = 1,3). Recall for any finite < U' < U, 
the system is always renormalized to the SU{4) symmet- 
ric model 19 ' 20) and the Kondo enhancement always appear 
at the gate voltage corresponding to Nd = 1,2,3; G(0, 0) is 
determined by the Friedel sum rule eq. (20) irrespective of 
the interaction asymmetry. We stress the interaction asymme- 
try manifests itself only at finite temperature and invisible at 
T = 0. It appears as a substantial 'dip structure' of G(T, 0) 
around half-filling. As we have claimed recently, 25 ' such fea- 
ture caused by interaction asymmetry reproduces quite well 
experimental observation in carbon nanotube dots 14 ' 15 - ) (see 
Sec. 7 for more details). 

To clarify how the asymmetric interaction affect the con- 
ductance profile particularly at half-filling, we make a direct 
comparison between the behaviors at Nd = 1 and Nd = 2 
in Fig. 3 by changing U'/U = 1, 0.9, 0.8, 0.6, and 0.4. By 
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Fig. 2. The temperature evolutions of the conductance G(T, 0) as a function 
of the gate voltage with varying asymmetry U'/U = 1, 0.8, and 0.6. 




Fig. 3. (color online) Temperature dependence of the conductance G(T, 0) 
as a function of the interaction asymmetry U'/U = 1, 0.9, 0.8, 0.6, and 0.4 
at Nj = 1 (quarter-filling) and at = 2 (half-filling). 



decreasing U' away from the symmetric point U' = U, the 
characteristic energy controlling thermal suppression seems 
reduced at Nd = 2, but enhanced at Nd= 1 ■ We will elucidate 
this tendency by estimating T* numerically in Sec. 4. 
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Fig. 4. Bias voltage evolution of the zero-temperature differential conduc- 
tance G(0, V) as a function of the gate voltage with varying the interaction 
asymmetry (a) U'/U = 1, (b) U'/U = 0.8, and (c) U'/U = 0.6. 



3.3 Effect of interaction asymmetry on G(0, V) 

We now turn our attention to the non-equilibrium transport, 
i.e., nonlinear conductance. To understand the interplay be- 
tween finite bias voltage effect and the interaction asymme- 
try, we here examine the zero-temperature differential con- 
ductance G(0, V). 

Figure 4 shows typical bias voltage evolutions of G(0, V) 
as a function of the gate voltage with asymmetry parameters 
(a) U'/U = 1, (b) 0.8, and (c) 0.6. We immediately notice that 
the profile of G(0, V) looks quite similar to that of G(T, 0). By 
applying larger bias voltage {eV/U = 0.03, 0.05, 0.1), the 
Kondo enhancement both at Nd = 1,3 and Nd = 2 is sup- 
pressed, but the latter suppression is more distinctive. Like 
the thermal evolutions, the interaction asymmetry manifests 
itself only at a fixed finite bias voltage, as a 'dip structure' of 
G(0, V) around half-filling. 

Figure 5 shows the zero bias peak G(0, V) either at N d = 1 
or at Nd = 2 with varying the asymmetry parameter U'/U = 
1.0, 0.9, 0.8, 0.6, and 0.4. Since we assume the symmetric 
coupling Ti = Yr to the leads, the G(0, V) becomes symmet- 
ric regarding V. From its behavior, we can read off how the 
characteristic energy scale depends on the asymmetric param- 
eter U'/U. The dependence is the same with what is observed 
in the temperature evolutions of G(T, 0); decreasing U' leads 
to the reduction of the characteristic energy at Nd = 2, but the 
enhancement at Nd = 1 . 




eV 



Fig. 5. (color online) The differential conductance G(0, V) as a function 
of the bias voltage with varying the interaction asymmetry U'/U = 
1,0.9,0.8,0.6,0.4 at N d = 2 (approaching G(0,V -» 0) -» 4G ), and 
at Nd = 1 (approaching G(0, V-t0)-) 2G ). 



Our numerical calculations show that the conductance at 
N d = 2 

4. Characteristic temperature T* 

The Kondo temperature is an energy scale characterizing 
the Kondo effect, i.e., the Anderson model eq. (1) at the cen- 
ter of an odd Coulomb blockade valley. It is a crossover tem- 
perature, not the transition temperature, so we can fix it only 
modulo numerical factor with ambiguity. Instead, we seek to 
define a characteristic energy scale at each gate voltage cov- 
ering the entire range of the topmost shell. 

Recalling the present effective theory expresses the dot 
Green function as 

G> = 0)= 1 =-2^), (21) 
~ Sd + iF ^s] ~ p 

we naturally define a characteristic temperature T* of the 
model 

r = Js 2 , + f 2 I (22) 

at each gate voltage or N d - We stress it is imperative to include 
the shift Sd in defining a consistent energy scale at each Nd. 

When Sd vanishes, we see T* = F, the Kondo peak width. 
This is a situation of the S U(2) Kondo effect. However, as 
derived from the Friedel sum rule of S U(4) model, eq. (20), 
the SU(4) Kondo peak needs to shift at the Nd = 1,3 17 - 18) 
by \s d \ = f = T*l V2. In contrast, the S U(4) Kondo peak at 
Nd -2 does not shift because of the electron-hole symmetry. 

The dependence of the characteristic energy T* obtained 
from our KR-SBMT calculations is shown in Fig. 6 for the 
entire range of the topmost shell (N d = ~ 4). Though the 
state degeneracy at U' — U is six-fold at Nd =2 and four-fold 
at Nd — 1 respectively, T* of the former is lower than that of 
the latter. 

The presence of the interaction asymmetry further ampli- 
fies this difference of T* as in Fig. 7; it reduces T* at Nd — 2 
but enhances it at Nd = 1 ; T* is almost unchanged at the filling 
corresponding to Coulomb blockade peaks (Nd ~ 1.3 and 2.7 
in Fig. 6). Because of this, the Kondo enhancement around 
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TAT* 



Fig. 6. Dependence of the characteristic temperature T* on the average 
electron number of the dot with varying interaction asymmetry V /U. 



Fig. 8. (color online) Universal scaling of the temperature scaled by T* 
Data are the same with Fig. 3. 




Fig. 7. (color online) Dependence of the characteristic temperature T* on 
the interaction asymmetry U'/U for Nj = 2 (filled square) and for Nj = 1 
(empty square). 



Nd » 2 is destroyed more rapidly than at Nd ~ 1 either by 
finite temperature or by finite bias voltage. The above depen- 
dence of T* gives a consistent explanation on the conductance 
profiles Figs. 2-5 in the previous section. 

How much asymmetry deforms the conductance profile 
of the S U{4) symmetry? The characteristic temperature T* 
at half-filling (of the symmetric model) controls it. When 
\U - U'\ exceeds T*(N d = 2, U' = U), it reduces T* rapidly, 
producing the conductance profile with a distinctive dip at 
half-filling. One can observe the effect only at finite temper- 
ature or by finite bias voltage. The above behavior of T* at 
Nd - 1,2 agrees qualitatively with other numerical estimates 
of the Kondo temperature. 19 ' 20 ' 35 ' 36 " 1 However, we need to dis- 
tinguish between the reduction of the characteristic energy 
and the crossover between the universality classes: the for- 
mer does not necessarily mean the latter. We will focus on the 
issue in next section. 

5. Universal scaling of linear and nonlinear conductance 

Having understood qualitative behavior of the conductance 
profiles by temperature or bias voltage, we now focus on 
examining an important issue: what is the universality con- 



trolling quantum transport in the presence of the interac- 
tion asymmetry? The nature at half-filling have been argued 
actively so far, proposing the S U(2) x S U(2) Kondo ef- 
fect; i9,20,35) or jjjg s u( 4) Kondo effect. 36) We should bear in 
mind that the absence of the exact S U (4) symmetry of the 
model does not necessarily legitimate the former universality, 
because the model flows toward the S U{4) symmetric strong 
coupling point for < U' < U; neither the reduction of T* in 
Fig. 7. 

When the nonlinear conductance G(T, V) is dominated by 
a unique universality, we expect to express G(T, V) by some 
universal function F„(t = T/T*,v = eV/T*) at each filling 
n = N d : 



G(T, V)/G = F n (t,v). 



(23) 



The Friedel sum rule eq. (20) constrains F n (Q, 0) = 
4 sin 2 (7rn/4). Otherwise, with a crossover present between the 
universality classes, we should observe a continuous change 
of the universal function F„{t, v). We test numerically the uni- 
versality by making universal scalings of the temperature evo- 
lutions and the bias voltage evolutions of the conductance. 

In the following, we will find the temperature or bias- 
voltage evolutions of the conductance G(T, V) indeed exhibits 
itself of the S U{4) behavior for different values of U'/U < 1. 
Hence we will claim the phenomena is dominated solely by 
the S U (4) universality, not a crossover between two different 
universality classes. 

5.1 Universal scaling of G(T,0) 

We first examine a universal form of the temperature evo- 
lutions of the linear conductance 



G(T,0)/G = F n (t,0); t = T/T*. 



(24) 



In Fig. 8, universal temperature dependence of G(T, 0) ei- 
ther at Nd = 1 or Nd - 2 are confirmed numerically, by 
rescaling Fig. 3 as a function of t — T/T* with varying 
U'/U = 0.4 ~ 1.0. Here we defined T* at each N d and 
U'/U. The results are striking. For all values of U'/U both 
at Nd = 1,2, universal curves collapse well up to t < 1 (an 
upper scale restricting the validity of the present analysis). It 
implies that the universal temperature dependence either at 
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Fig. 9. (color online) Universal scaling of the bias voltage scaled by T" 
Data are the same with Fig. 5. 



half-filling or at quarter-filling reduces to the S U (4) symmet- 
ric case. Additionally Fig. 8 shows that such universal depen- 
dence differs slightly but significantly between Nj = 1,2. The 
universal curves obtained here look very similar to those of 
the symmetric model U' — U obtained in NRG. 32) 

5.2 Universal scaling o/G(0, V) 

We now turn our attention to universal scaling of non- 
equilibrium transport regime, of the nonlinear conductance 



G(0, V)/G = F„(Q, v); v = eV/T* 



(25) 



In Fig. 9, we show the bias-voltage evolution G(0, V) as a 
function of the scaled bias voltage v = eV/T* either at = 
1 or at Nd = 2. As the most important observation, we can 
confirm that all of the universal scaling curves for U'/U = 
0.4 ~ 0.9 collapse on the S U (4) symmetric one, similarly to 
temperature dependence. In comparing with thermal universal 
curves F„(t, 0) in Fig. 8, bias-voltage universal curves F n (Q, v) 

at Nd = 1,2 exhibit rather similar behavior, up to t < 1 (T < 

n 

5.3 Universal scaling of G(T,V) 

As a concrete dependence, thermal universal curve F„{t, 0) 
and bias-voltage one F „(0, v) behave differently for n = Nd — 
1,2. This means that we need three independent variables f, 
v, and n fully to describe a universal curve. However, a sim- 
pler description is possible to characterize the suppression by 
finite temperature or bias-voltage. The idea essentially comes 
from seeking what is a cut-off energy of the RG flow. 

A typical nonlinear conductance evolution at finite tem- 
perature is shown in Fig. 10 for U'/U = 0.6, for (a) T = 
0.0025 U, (b) 0.025 U and (c) 0.05 U. In this system, the 
characteristic temperature T* at half-filling is estimated as 
T*(Nd = 2) = 0.023t/. We understand its behaviors by classi- 
fying these behaviors, comparing with this energy scale T*(2). 
At (a) T = 0.0025 U <sc T*(2) (t <sc 1), the conductance 
profile looks similar to that at T = in Fig. 4 (c). It shows 
the temperature is so small that no thermal effect appears. 
On the other hand, at (c) T = 0.051/ » T*{2) (t » 1), 
all the conductance profiles with different bias-voltage V al- 
most collapse except for eV = Q.1U which is larger than the 
temperature T = 0.05 U, which is a hallmark of the temper- 
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Fig. 10. Conductance profiles G(T, V) by varying different temperature and 
bias voltage, (a) T = 0.0025(7 (b) T = 0.025(7 (c) T = 0.05(7. For each 
temperature, results of bias voltage eV = 0.005, 0.03, 0.05, and 0.1(7 are 
shown. The characteristic energy at half-filling is 7**(2) = 0.023(7. 



ature T » T*(2) controlling the conductance profile for a 
region T > eV. Consistent behaviors are also confirmed at 
(b) T = 0.025 U ~ T*(2): the temperature controls the con- 
ductance profile for T > eV; the bias voltage, for T < eV. In 
Fig. 11, we illustrate universal functions F n (t,v) for n = 1,2 
of the S (7(4) universality, obtained by our KR-SBMT ap- 
proach. We conclude that the larger energy scale max[f, v] 
between the temperature and the bias voltage effectively dom- 
inates the nonlinear conductance profile. The result presented 
in Fig. 11 is indeed calculated for U'/U = 0.6, but as we 
have shown so far, this universal function F„(t, v) is almost the 
same up to T* with that of U'/U = 1 and also of U'/U < 1. 

Before concluding this section, we give two cautionary 
comments regarding the applicability of the present univer- 
sal dependence of the conductance: (1) on conductance pro- 
files observed as a function of the gate voltage e^; (2) on the 
dependence on the value of U/T. The first point is related to 
our employed approximation (KR-SBMT). Although we be- 
lieve the present analysis of the T* scaling may be useful even 
for max[T, e V] > T*, our present KR-SBMT method is be- 
lieved to describe the system adequately only up to < T*; it 
does not fully incorporate charge fluctuations, which will be 
important for transport under a large bias. The characteris- 
tic temperature T* is lower at half-filing, higher at quarter- 
filling, particularly in the presence of the interaction asym- 
metry (Figs. 6, 7). Hence the conductance evaluated by the 
KR-SBMT is validated up to a smaller value T or eV at half- 
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Fig. 1 1 . (color online) Universal function F„ ( t, v) of the S U (4) universality 
class, obtained by the KR-SBMT. (a) G(7\ V)/G(0, 0) = F(t, v)/4 at half- 
filling (N d = 2) and (b) G(T, V)/G(0, 0) = Fit, v)/2 at quarter-filling (N d = 
1), 



filling than at quarter- filling. Such signature is observed in 
the conductance profile for eV = Q.1U in Figs. 4 and 10, 
where max[T, eV]] » T*(Nd = 2) and the Coulomb blockade 
physics dominates the system. In those cases, the conductance 
is anticipated to vanish, with some nonlinear bias effect that 
the KR-SBMT is likely to miss. The second point is related to 
why our analysis observes the S U (4) universality class, disre- 
garding the interaction asymmetry, without a crossover. This 
important issue will be examined further in next section. 

6. T* renormalization or a crossover: Effect of finite U 

As we have shown, our results of conductance at finite 
temperature and/or bias voltage behave consistently with the 
S U{4) universality class even in the presence of the interac- 
tion asymmetry U'/U < 1, exhibiting no crossover toward a 
different universality class. The results may seem odd at first 
glance, contradicting with the suggested crossover between 
SU{4) and SU(2) x SU(2) universality classes. 19) However, 
such a crossover behavior takes place only for a larger value 
of U/T. In a usual range of experimentally feasible parame- 
ters, we claim the renormalization of the characteristic tem- 
perature T*, not a crossover of the universality classes, plays 
a more dominant role, controlling the conductance behavior. 

To support our view, we also examine the behavior with 
larger values of U. Figure 12 shows the characteristic temper- 
ature T* at Nd — 2 as a function of the interaction asymme- 
try U'/U for U/T = 50,60, to be compared with Fig. 7. In 
this case, the interaction asymmetry U'/U reduces T* quite 
rapidly, which signals a crossover of the universality class 
from SU(4) to SU{2) x SU(2). l9) We confirm that such a 
crossover indeed occurs by performing the universal scaling 
(Fig. 13). For U/T = 50, the universal scaling indicates that 
the scaled curve for U' /U — 0.9 is clearly different from that 
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Fig. 12. Dependence of the characteristic temperature T* on the interaction 
asymmetry U'/U at N d - 2 for U/T = 50 and 60. 
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Fig. 13. Universal scaling of the bias voltage scaled by T* at Nj = 2 for 
U/T = 50. A crossover between different universality classes occurs. 

for U'/U = 1 or 0.99; the universality class crosses over be- 
tween S 1/(4) and S U{2) x S U(2) in this case. 

In usual experiments, however, U is roughly in the range 
of U = 10 ~ 20F (see also the next section). For this case, 
as was shown in §5, the dominant effect at Nd ~ 2 is due to 
the renormalization of T*, not a crossover of the universality 
classes. 

7. Experimental manifestation 

We have found theoretically the effect of the interaction 
asymmetry on the conductance profile substantial, either at 
finite temperature or at finite bias voltage. Currently, a few ex- 
perimental realization of the S t/(4) Anderson model eq. (1) 
are known: transport through a carbon nanotube dot, and, 
through double quantum dots coupling capacitively with each 
other. 3, 4) We discuss our present results in the light of existing 
experimental data. 

7.1 Carbon nanotube dot 

Several groups have investigated and observed the Kondo 
enhancement in the linear conductance measurement through 
a carbon nanotube dot. 12 ~ 15) A carbon nanotube dot has dou- 
bly degenerate orbitals in the topmost shell, and moreover, 
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Fig. 14. (color online) Schematic plotting for the comparison with carbon 
nanotube dot experiments: Conductance as a function of the gate voltage 
for U'/U = 0.997 and U = 30r. \U - U'\ « T*(N d = 2, U' = U) is 
realized. Inset shows larger coupling parameter regime U = 15T. 




Fig. 15. (color online) Schematic plotting for the comparison with carbon 
nanotube dot experiments: Conductance as a function of the gate voltage 
for U'/U = 0.9. Other parameters are the same with Fig. 14. \U - U'\ » 
T*(N d = 2, V = U) is realized. 

the conductance according to the S U (4) Friedel sum rule 
eq. (20) is clearly observed experimentally at low tempera- 
tures. 15 ' This observation clearly justifies the validity of the 
S t/(4) Anderson model. Notwithstanding, we have not under- 
stood fully the temperature evolution of the conductance, par- 
ticularly around Nd ~ 2: While the experiment by Makarovski 
et al. has observed visible low-temperature "Kondo-like en- 
hancement" around Nd = 2, 15) such enhancement was absent 
exhibiting a characteristic 'dip' in others. 13 ' 14 - 1 We suggest 
the interaction asymmetry enables us to understand system- 
atically finite temperature evolutions in these experiments. 25) 
In Figs. 14 and 15, we present schematic calculations mim- 
icking the experimental situation, changing U'/U slightly, 
with U = 30r (with U = 15r for the insets). In Fig. 14 with 
U'/U = 0.997, the conductance profile reproduces all the 
features of the S U (4) symmetric model. Notably, the Kondo 
enhancement at low temperature occurs either at Nd ~ 2 
or at Nd ~ 1 similarly. As was claimed already in ref. 32, 
these temperature evolutions by the S U(4) symmetric Ander- 
son model agree very well with what is observed either at 
V g ~ 3.9V or at V g ~ 5.3V in the experiments by Makarovski 
et al. I5) To our surprise, however, the conductance profiles 



with four peaks (U - 30T) are modified considerably by a 
relatively small asymmetry (less than 10%), having a charac- 
teristic dip structure around half-filling as in Fig. 15. This is 
because the characteristic energy at Nd -2 is rather small, 
T*(N d = 2,U' = U) * 0.014C7. The value U'/U = 0.9 
conforms to almost equally-spaced Coulomb blockade peaks, 
observed at T — 8K. I4) We believe Fig. 15 captures essential 
characteristics experiments 13 ' ' with a reasonable choice of 
physical parameters; for instance, in Fig. 3(a) or Fig. SI2 of 
Jarillo-Herreo et al., ^ the conductance at V g ~ 3V is smaller 
than that of V g ~ 2.8V or 3.2V at each temperature. 

Regarding the mixed-valence regime {U/T = 15, insets of 
Figs. 14 and 15), four peaks merge to form one big peak show- 
ing the Kondo enhancement. In this regime, we find the in- 
teraction asymmetry U'/U hardly affecting the conductance 
profile even at finite temperature. 

7.2 Double dot system with capacitive coupling 

Several groups have been investigating experimentally on 
orbital and/or spin Kondo effects in double quantum dot 
systems with capacitive coupling, 3 ' 4 ' 7 ' 8) where the orbital 
Kondo effect has been observed at quarter-filling. The mea- 
sured Kondo temperature is found surprisingly high. Thus 
we should possibly attribute the phenomena to the SU(4) 
Kondo effect, as an entanglement of spins and orbits. There- 
fore, within the universality of the topmost shell, we will be 
able to describe the system by applying the S U (4) Anderson 
model eq. (1). 

Although the S U(4) Anderson model suggests that the 
conductance should exhibit the Kondo enhancement both at 
quarter-filling and at half-filling, to the best of our knowledge, 
no enhancement has been observed so far at half-filling. How- 
ever, we argue that the absence of the Kondo-like enhance- 
ment at half-filling is due to the combined effect of finite bias 
effect and the interaction asymmetry in those systems. Con- 
sequently, it is expected that the conductance enhancement 
should appear also at half-filling by suppressing the bias volt- 
age at low temperature. 

To make our argument more concrete, let us examine typ- 
ical double dot experiments, 3 4) whose IV characteristics at 
T = 30mK are excerpted in Fig. 16 (a), which is mea- 
sured by applying the source-drain voltage 80/iV. The volt- 
age V(i-2) controls the degeneracy of single levels; the top- 
most single orbits realize two-fold degeneracy at V(i-2) ~ 
3.75,4.25,4.75,5.25 etc. The honeycomb structure clearly 
shows Kondo enhancement at quarter-filling (denoting re- 
gions 'b' in Fig. 16(a)); no enhancement at half-filling. In 
Fig. 16 (b), choosing the parameters U = 0.7meV, U'/U - 
0.6, and U/T = 20, we plot corresponding KR-SBMT results 
of the nonlinear conductance G(T, V) by varying the source- 
drain voltage V = 0,35,70yuV. The result of V = 70yuV re- 
produces well with experimental data; Kondo enhancement is 
observed at quarter-filling but not at half-filling. 

Our KR-SBMT calculations evaluate the characteristic en- 
ergy at half-filling T*(2) to be 16/^eV for these systems. So 
we conclude experiments realize large bias regime suppress- 
ing the Kondo effect at Nd ~ 2. It also suggests that when 
we can decrease the source-drain voltage, less than 16//V, the 
Kondo-like enhancement should appear at half-filling as well 
as at quarter-filling. 
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Fig. 16. (a) Experimental IV characteristics excerpted from refs. 3,4 (see 
the text), (b) A schematic differential conductance obtained by the KR- 
SBMT at T = 30 mK, as a function of the gate voltage, varying bias voltage 
V = 0, 35, 70/iV. 



8. Conclusion 

We have investigated the role of the inter- and intra-orbital 
interactions in transport through quantum dots with two-fold 
orbital degeneracy, e.g., carbon nanotube dots, double dot 
systems with capacitive interaction etc. By using the SU(4) 
Anderson model with the interaction asymmetry between or- 
bits/dots and within the validity of the KR-SBMT, we have 
shown how significantly a small amount of it can affect the 
profile of the linear conductance at finite temperature, and of 
the nonlinear conductance. We have shown that the phenom- 
ena can be understood systematically by the different depen- 
dence T* on interaction asymmetry at each gate voltage. In 
addition, we have compared our theoretical results with exper- 
imental data in two typical situations: carbon nanotube dots, 
and double dot systems with capacitive interaction. We be- 
lieve our theoretical results agree very well, suggesting unex- 
plored phenomena of the Kondo enhancement in double dot 
systems. We have shown clearly the interaction asymmetry 
is essential in explaining conductance profiles through these 
two-fold degenerate systems in real experiments. 
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